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Abstract 



We compute the threshold-resummed cross section for pseudo-scalar MSSM Higgs boson production by gluon 
fusion at hadron colliders. The calculation is performed at next-to-next-to leading logarithmic accuracy. We 
present results for both the LHC and Tevatron Run II. We analyze the factorization and renormalization scale 
dependence of the results, finding that after performing the resummation the corresponding cross section can be 

computed with an accuracy better than 10%. 
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1 Introduction 



The hunt for the Standard Model (SM) Higgs boson is clearly one of the biggest physics goals at the LHC. Its 
search is not only a great challenge from the experimental point of view, it also requires a huge theoretical effort 
in order to provide very precise predictions both for signal and background. Besides the SM, there are many 
other options for New Physics (NP). Among the possible scenarios, the MSSM is a promising one. It provides a 
solution for the hierarchy problem and also introduces a good dark matter candidate, the lightest supersymmetric 
particle. The Higgs boson sector of this model consists in two complex Higgs doublets. After the EWSB, three 
degrees of freedom are absorbed by the electroweak massive bosons and the remaining five give rise to the usual 
SM scalar Higgs boson (/i), a heavier neutral one {H), two charged scalars (H^) and a pseudoscalar neutral Higgs 
(A). In the MSSM, the tree level masses depend upon two parameters, which can be selected to be the mass of the 
pseudoscalar Higgs (m^) and the ratio of the vacuum expectation values of the two doublets tan f}= V2/V1 [1,2]. 
These parameters are constrained by LEP and Tevatron experiments, with bounds given by tua > 100 GcV and 
with the interval 0.5 < tan (3 < 2.4 currently excluded [3]. The actual status of SM and SUSY Higgs production at 
LHC is thoroughly reviewed in Ref. [4]. 

The main channel for the production of an on-shcU neutral Higgs boson at the LHC, both for SM and MSSM 
bosons, is the gluon fusion process. Therefore, the corresponding cross section needs to be under good theoretical 
control. The calculation, for both CP-even and CP-odd Higgs particles, has been performed up to ncxt-to-lcading 
order (NLO) in an exact way [5]. Surprisingly, the NLO corrections turned out to be quite large, accounting for 
an increase of almost 100 % of the leading-order (LO) result. Furthermore, adding the NLO term did not result in 
a sizeable reduction of the renormalization and factorization scale dependence. Nevertheless, the NLO corrections 
showed that the gluon-Higgs interactions could be very well approximated by using an effective Lagrangian ap- 
proach, inspired by low energy theorems. This effective theory is obtained by considering an infinite mass limit for 
the top quark [6]. This approximation allowed to perform the next-to-next-to leading order (NNLO) calculations 
for the inclusive cross section and even to analyze the transverse momentum distribution of the Higgs boson [7,8]. 

At NNLO accuracy, the dominant soft- virtual corrections were first computed in [9,10] , while the full result [11- 
14] became available later. For pseudoscalar Higgs, the full NNLO contribution was computed in Refs. [11,15,16]. 
These corrections turned out to be considerably smaller than the NLO ones, bringing a reasonable reduction in the 

scale dependence as well. 

The soft-gluon resummation, which includes in the cross section the dominant effects of the higher order terms 
of the perturbative QCD series, was performed only for the scalar Higgs, up to next-to-next-to- leading logarithmic 
(NNLL) accuracy [17]. These corrections showed a modest numerical impact, and a further (slight) reduction of the 
scale dependence, providing a hint on the convergence of the asymptotic QCD perturbative series. The quantitative 
reliability of the soft-gluon approximation was probed by comparing the truncation of the resummed result at NLO 
and NNLO with the fixed-order calculation [10,17]. This comparison showed that the resummed result can also 
be trusted away from the threshold region, thus convalidating the accuracy of the threshold resummation in the 
entire kinematical phase space. 

In the present work we compute the cross section for CP-odd Higgs production by gluon fusion up to NNLL, 
in the infinite top mass limit. As the corrections for scalar and pseudoscalar Higgs do not greatly diff'er, we still 

expect, based on the scalar Higgs result, that the soft-gluon resummation is a highly accurate approximation to 
the exact result in the whole phase space, and not only around the threshold region. 

This article is organized as follows. In Section 2 we introduce our notation and briefiy review the theoretical 
framework, including QCD cross sections at fixed-order and the basics of soft-gluon resummation. In Section 3 
we present and discuss our results, both for the LHC and the Tevatron Run II while in Section 4 we give our 
conclusions. 
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2 Theoretical frame 



According to the mass factorization theorem, the inclusive cross section for the production of the pseudoscalar 
Higgs boson by the collision of hadrons hi and ft,2 may be written as 



Cr(s,M|)=V'/ dxi dx2 fa/hi{xi,fJ,F) fb/h2ix2,fJ-%) dz < 

„ u Jo Jo 



TA 

XIX2 



a,b 

■aab{s,Ml){z;as{fil),Ml/f,l;Ml/fil) , (1) 

where Ma is the pseudoscalar Higgs boson mass, ta — Al'^/s, and and fin are the factorization and renormal- 
ization scales, respectively, dab denote the partonic cross section for the process a + b ^ A + X , computable in 
perturbative QCD. The parton densities of the colliding hadrons are denoted by /a/?i(a^7 /^/p), where the subscript a 
labels the parton type. We use parton densities as defined in the MS factorization scheme. For practical purposes, 
one works with the hard coefficient functions Gab instead of the partonic cross sections, the first one given by 



fafc 



a^'^zGabiz), (2) 



where (t(0) is the Born level contribution. The incoming massless partons a, b do not couple to the pseudoscalar 
Higgs boson directly. In hadron collisions the main production mechanism is through heavy quark triangle loops 
and therefore, the total cross section also depends on the top (Afj) and bottom {Mi,) quark masses. 

The NLO coefficients G^^^^ have been exactly computed in Ref. [5], where it was also observed that the NLO Higgs 
boson cross section can be well approximated in the low tan f3 regime by considering its limit Mt 3> Ma [18]. Hence, 
along this paper we will work within the large- Afj approximation: we consider the case of a single heavy quark 
(the top), and Nf = 5 light-quark flavors, neglecting all the contributions to G^^^^ that vanish when Ma/Mi 0. 
Nevertheless, the full dependence on Mt and M^, is included in cr*^°^ in order to improve the accuracy of the 
calculation. The large- Mj approximation allows the use of the effective-Lagrangian approach [6,19,20], that 
shrinks the top quark triangle loop into an effective point-like vertex, considerably simplifying the evaluation of 
the Feynman diagrams involved in the process. 

The heavy top mass limit is reliable as long as ta,v? f3 <C Mt/Mb [21], because of the quark-Higgs couphngs. 
Therefore, any calculation relying in the Mt 00 approximation is valid for low values of tan /3 . Moreover, 
we assume that the lighter squarks are much heavier than the top quark, thus their contribution can be safely 
neglected. 

The coefficient function Gabiz) is dominated by soft terms in the limit z ^ 1. Therefore, our main objective is 
to study the effect of soft-gluon contributions to all perturbative orders. This task requires to work in the Mellin 
(or A^-moment) space [22,23]. We thus introduce our notation in the iV-space. 

We consider the Mellin transform aN{Mj^) of the hadronic cross section cr(s, Mj^). The A^-moments with respect 
to Ta = M\/s at fixed Ma are customarily defined as follows: 

(Tm{MI)= I dTAT^-' a{s,Ml) . (3) 
Jo 

In A^-moment space, Eq. ([1]) takes a simple factorized form 

(Tn{MI) = fa/huNilJ^-) fb/h2,N{fJ-%) CFab..N{.as{^i]i),M\/ ii%;M\/ ^i],) , (4) 

a.b 

where fi/h,N a-nd dij.N represent the Mellin transforms of the parton distributions fi/h and of the partonic cross 
sections crij, respectively. 

In order to perform the threshold resummation, we first note that the threshold region z — > 1 corresponds 
to the limit A — > cx) in Mellin space. The dominant contribution in this limit is due to the large logarithmic 
terms Og In"* A . Being the only channel open at the Born level, the gg ^ A subprocess is the unique partonic 
contribution that can give rise to the large threshold logarithms. The formalism to systematically perform soft- 
gluon resummation for hadronic processes was developed in Refs. [22,23]. In the case of Higgs boson production. 
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one has 

[' +c>o 2n 



G,,, ^ = 4 1 + E "5 E Gr""^ In" iV + 0(1/7V) = g(';1 + 0{1/N) , (5) 



n— 1 m— ^ 

where the dominant contributions in the large- iV limit may be reorganized in the following all- order resummation 
formula: 

■ eMGhiasi^^l)^nN■,Ml/^ll,Ml/^^l)} . (6) 

The function Cgg{as) contains all the contributions that are constant in the large- iV limit. They are originated 
from the hard virtual contributions and non-logarithmic soft corrections, and can be computed as a power series 
expansion in as- The large logarithmic terms In™ TV (with 1 < m < 2n), which are due to soft-gluon radiation, 
are included in the exponential factor expC//i. It can be expanded as 

gn (as{til)MN- Ma Mh_\ InTV .9«(A) + .gf (A, MIK; M1/^^) 

+ o.s{A) (A,Ml^?j;Ml^^) -f 0{as{as\nNf) (7) 
where A — bQas{fJ.%) InN and bo is the first coefficient of the QCD /3-function. 

Due to the universality of the soft-gluon emission, the Q factor is independent on the type of Higgs boson 
produced in the final state. Hence, the coefficients of the expansion are the same for both h and A. The expressions 
for these coefficients can be found in Ref. [17]. The difference in the resummed expansion between the scalar 
and pseudoscalar Higgs only shows up in the Cgg factor which, being partially originated on the hard-virtual 
contribution, obviously depends on the Higgs type under consideration. For the sake of brevity, we only write the 
differences between the A and the h terms, which read 



where Nf is the number of different light quark flavors. 

Going back to Eq. ([7]), the term InA^ g^p resums all the leading logarithmic (LL) contributions In""*"^ A^, 

(2) (3) 

gj^ contains the next-to-leading logarithmic (NLL) terms Og In" N , asgh collects the next-to-next-to-leading 
logarithmic (NNLL) terms a^^^ In" N , and so forth. In this context, the product 015 IniV is formally considered as 
being of order unity. Therefore, the ratio of two successive terms in the expansion ([7]) is formally of 0{as)^ which 
makes the resummed logarithmic expansion in Eq. ([7]) as systematic as the usual fixed-order expansion in powers 
of as. 

The leading coUinear contributions can also be included in the soft-gluon resummation formula by performing 
[17, 19] the following shift 

C'^J^CfJ^^'^ (9) 

that correctly resums all the terms of the type (a^ In^"^^ N)/N that appear in Ggg"* . 

When attempting for the resummation, one is interested in taking some advantage of the full fixed-order cross 
section calculation as well. It is therefore customary [17] to perform a matching between both approaches, which 
can be schematically written as 

^matched ^res i ^f.o ^res\ { A r\\ 

<T =(7 + a' ~ a l/.o (lUj 

where a^'^'^ corresponds to the result obtained using Eq. ([6]), a^-° is the fixed-order cross section and cr'"'^'*|/.o 
represents the expansion of the resummed result at the same order in ag as the fixed-order result. This improved 
matched cross section is our final result for the process. In consequence, throughout the paper we shall refer to 
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the different orders of the matched cross sections directly as LL, NLL and NNLL. The accuracy of the matching is 
assured by the order at which the Cgg coefficient is computed, being obtained by a direct comparison between the 
resummed and the fixed-order calculation. Therefore, in a^ " — cr''^*|/.o only the hard terms, which are strongly 
suppressed in the N oo limit, survive. 

Recently, the function ^jj*' was presented [24,25]. In principle, it allows to perform the resummation up to 
NNNLL accuracy. However, the full matched calculation at NNNLL, can only be done if the fixed-order NNNLO 
result [3 were available. Nowadays, only the soft contribution was derived [25]. We have decided to look at the 
effect of including the g^^ term and setting Cg^^ = 0. This inclusion leads to a really slight modification of the full 
NNLL result thus validating the convergence of the resummed series expansion, and allowing us to safely neglect 
the g^'^ function along this work. 

Finally, once the expression of the cross section has been computed in iV-space, the physical result can be 
obtained by Mellin inversion. In order to avoid the Landau singularities explicitly present in the exponential factor 
in Eq. ([7|) we use the Minimal Prescription, as described in Ref. [26]. 



3 Results 

We have developed the program THIGRES, a FORTRAN code to compute the resummed (fixed-order) cross 
section up to NNLL (NNLO) accuracy, both for scalar and pseudoscalar Higgs boson in the heavy Alt limit. The 
improvement over previous calculations lies in the fact that the partonic cross sections up to NNLO are directly 
written in Mellin space. The Mellin coefficients of the hard functions G^^'^^ were presented in [27]. We have 
recomputed these coefficients in an analytical way, by performing the Mellin transform of the results presented in 
Refs [11-13,15,16], finding some non-negligible differences in the gg and qq channels at NNLO with respect to [27]. 
The essential ingredients for Mellin transformation can be found in [28] . We have taken advantage of the ANCONT 
Fortran code [29] which provides most of the required special functions in A^-space. 

One essential missing element to tackle the calculation in Mellin space are the PDFs. The available parton 
distributions are always given in the x space. In order to transform them to A'^-space, we first perform a fit of 
the densities at the needed scale using a functional form that allows for a simple analytical Mellin transform. We 
find that a linear combination of Eulerian functions x^"\l — x)^ , which in Mellin space give rise to beta functions 
B{a + N, f3 — 1) are enough to reproduce all the features of the usual parton distributions [30]. After performing 
some clever sampling, the result of the fit allows to compute A^-moments analytically with an accuracy better than 
0.5%. Having both the necessary fixed-order and resummed coefficients and the PDFs in Mellin space, one is able 
to perform the calculation in the most efficient way with a considerable reduction in the required computer time 
and a gain in precision. 

We have worked with the MRST set of PDFs; using the 2001 LO [32] and the 2002 NLO and NNLO [33], 
although the code THIGRES allows the use of other PDFs. For the presentation of our results we use Mt — 176 
GeV and Mb = 4.75 GeV. Therefore, the cross section for the pseudoscalar Higgs boson is reliable only if tan /3 < (jl. 

In order to check the validity range the approximation in tan /? , we will start studying the dependence of 
the results upon this parameter. As we are neglecting the finite top mass effects in the partonic cross sections, 
tan /3only appears in the Born term (Tq. In Fig. [l]we plot ctq as a function of tan (3 , showing the corresponding 
variation of the Born cross section, for a AIa = 115 GeV boson at the LHC, according to whether one includes 
top and/or bottom mass effects. It is clear that the use of the infinite top mass limit in the Born cross section 
is not reliable. At least the Mt dependence must be included in ao- The solid (top -I- bottom) and dashed (only 
top) curves are similar around tan 13= 1; for this particular value, the pseudoscalar Higgs coupling to up and 
down type quarks is the same as for the scalar one. The inclusion of the bottom mass can be neglected until one 
enters the region where tan /3 becomes close to ^/ Mt /M^ Near this point, both curves start to separate from 
each other. It is rather noticeable that for values of tan /3> 10 the bottom effects are completely dominant, and 
therefore any calculation relying on the infinite top mass approximation can not be trusted for those values. This 

tAt least, the full soft- virtual contributions are necessary to compute the coefficient Cgg' . 

^We have explicitly checked the accuracy of the infinite top mass approximation by comparing our NLO results with the exact NLO 
ones provided by the FORTRAN code HIGLU [31]. Within this bounds the accuracy is always better than 10 percent. 
5 At this particular value, the A-top coupling is equal to the A-bottom coupling. 
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Figure 1: tan P dependence of the Higgs production Born cross section at the LHC for ~ 115 GeV, in the 
infinite top mass limit (dots), including top mass effects (dashes) and both top and bottom mass effects (solid) . 



behaviour is certainly expected, since the coupHngs of the pseudoscalar Higgs boson to the up (down) type quarks 
are suppressed (enhanced) by a factor tan (3 . 

It is very important to know the dependence of the cross section upon the renormahzation and factorization 
scales, as a way to estimate the size of the higher order terms not yet included in the perturbative expansions and 
therefore evaluate the uncertainties on the theoretical calculations. This dependence is shown in Fig. [51 considering 
a Ma = 150 GeV Higgs boson at the LHC, using tan /3= 3, for the fixed-order LO, NLO and NNLO results. In 
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Figure 2: Scale dependence of the Higgs production cross section at the LHC for Ma ~ 150 CeY, tan /?= 3, at LO 
(dots), NLO (dashes) and NNLO (solid). 

this figure, both scales were varied from Ma/4: to 4 A/^ in three different ways. In the plot on the left, the varying 
scales were chosen to be equal {/j-f — = X^'^a)- In the plot on the center the factorization scale was changed 
and the renormahzation scale was kept fixed {hr = Ma^^j^f — X^a), while in the one on the right, the opposite 
variation was performed (/ii? = Ma, piR = x^a)- As expected from the running of as, the cross section typically 
decreases when ^r increases. This effect is clearly noticeable in the right-side plot, and, more moderate in the 
left-side plot. The graph on the center shows how the variation oi hf leads to and opposite behavior, ie: the cross 
section increases with the growing of /ij?. This can be explained by the following fact: at the LHC the cross section 
is mainly sensitive to partons with momentum fraction x ~ 10~^. In this x range, the scaling violation of the parton 
densities is moderately positive and therefore one observes an artificially reduced factorization scale dependence. 
In the left-side plot one sees a partial compensation of the two effects, although the piR variation clearly dominates. 
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Figure 3: Scale dependence of the Higgs production cross section at the LHC for Ma — 150 GeV, tan P— 3, at LL 
(dots), NLL (dashes) and NNLL (solid). 



Another interesting feature of this plot is that it shows how sizeable are the higher order corrections. The change 
when going from LO to NLO is quite large, while the inclusion of the NNLO corrections has a moderate impact. 
We can consider this fact as a hint for the convergence of the perturbative series. Fig. [3] shows the same as Fig. [2l 
for the resummed cross section. One sees that the plot on the right keeps the typical dependence on fin, due to the 
running of the coupling constant. Nevertheless, the behavior for fixed fijj has changed, especially if we compare the 
higher orders (NLO vs NLL and NNLO vs NNLL). The rather flat result of Fig.[2]is now replaced by a considerably 
higher variation, due to the inclusion of both soft and collinear higher order terms. Therefore, in the left-side plot, 
the scale variations are fairly compensated, and particularly the NNLL result exhibits a tenuous dependence on 
the combined scales. Fig. [5] shows that the scale dependence is very slightly reduced when going from LO to NLO 
(as was already mentioned), and considerably reduced when going from NLO to NNLO. The implementation of 
resummation effects leads to a further reduction of the scale dependence. 

In Fig. [4]we present the fixed-order scale dependence for Ma = 150 GeV with tan (3=5 now at the Tevatron. 
In this Figure, many of the overall features that appeared in the LHC plots are present. The right plot still shows 
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Figure 4: Scale dependence of the Higgs production cross section at the Tevatron for Ma=150 GeV, tan /3= 5, at 
LO (dots), NLO (dashes) and NNLO (solid). 

the same dependence, again because of the running of ag. In the middle plot, in direct contrast with the result 
from Fig. [21 the cross sections increases with fip. At the Tevatron, the partons with roughly x ^ 0.1 contribute 
the most and in this region the scaling violation is slightly negative. Then, in the first plot both scales lead to an 
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overall decrease in the cross section. Fig. [S] shows now the scale dependence of the resummed cross section with 
the same parameters of the previous figure. The results are very similar to those in Fig. [31 Here, it is noticeable 
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Figure 5: 5*00^6 dependence of the Higgs production cross section at the Tevatron for Ma = 150 GeV, tan /?= 5, 
at LL (dots), NLL (dashes) and NNLL (solid). 



the rather flat scale dependence of the NNLL result in the left plot. It is quite remarkable that the LO and LL 
curves look very much alike. The effects of threshold resummation become important only when going to higher 
(NLL and NNLL) orders. As an overall feature of the scale dependence graphs, it is important to stress the fact 
that the resummed cross section is clearly more stable against scale variations than the fixed-order result. 

The importance of higher-orders effects in commonly presented through the introduction of K-factors, which are 
defined as the ratio of the cross section at a given order over the LO result. As it was mentioned before, within the 
large- Mt approximation, the higher order cross sections are proportional to uo, which is the only term that depends 
upon tan (3 . Therefore, in the infinite top mass limit, the K-factors are fully independent on that parameter. In 
Fig. [S] we present the K factor at LHC including its scale dependence. The bands were obtained by independently 
varying the scales in the region 0.5M^ < fj,p, < 2M^, with the constraint 0.5 < /iF/^fl < 2. The LO result that 
renormalizes the K factor was computed with fip = fiR = Ma- The NLO K factor is around 2, accounting for an 




Figure 6: Fixed-order and resummed K factors for Higgs production at LHC. 



increase in the cross sections of approximately the same amount as the LO result itself. The NNLO K factor shows 
a rather more moderate impact. The inclusion of soft-gluon effects at NLL and NNLL accuracy slightly increases 
the cross section on top of the fixed-order contributions and show a larger overlap between the corresponding bands, 
indicating a better convergence for the resummed series. It becomes clear the reduction of the scale dependence 
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for the higher orders, as the bands become thinner as the order grows. We also notice an increase of the K factors 
with Ma, consistent with the fact that the soft-gluon contributions become more dominant as the process gets 
closer to the hadronic threshold. Once the resummation is performed, the uncertainty due to scale variation can 
be estimated to the order of 10 percent. 

Finally, Fig. [7] shows the K-factor Higgs mass dependence at the Tevatron. Here we can see the same overall 
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Figure 7: Fixed-order and resummed K factors for Higgs production at Tevatron. 

features as presented in Fig. [6l As a major difference, we note that the K factors are considerably bigger for the 
Tevatron. This is due to the fact that the Tevatron center of mass energy is closer to the hadronic threshold, which 
is the kinematical region where soft-gluon effects become relevant. 



4 Conclusions 



We have presented the resummed cross section for pseudoscalar Higgs production by gluon fusion at NNLL accuracy 
in hadronic colliders, presenting the most relevant results for the LHC and Tevatron. 

Comparing with the fixed-order calculation, the numerical impact of the resummation effects was found to be 
rather moderate. This slight variation when going from NNLO to NNLL is providing a hint for the (hopefully) 
faster convergence of the perturbative series. Probably the most striking feature of the resummed result is the 
considerably reduction of the scale dependence. It allows to make theoretical predictions with a precision of about 
10 %, which are accurate enough for discovery of a Higgs boson at the LHC and Tevatron. Moreover, our results 
presents a probe to supersymmetry, as can also be useful to directly test the MSSM and/or another supersymmetric 
extension of the Standard Model in the low tan [3 regime. 

The Fortran code THIGRES, which computes total cross sections for both scalar and pseudoscalar Higgs boson 
up to NNLO (or NNLL), is provided upon request from the authors. . 
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